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ABSTRACT 

(D 

' We follow the timing properties of the neutron star low-mass X-ray binary 

system 4U 1705-44 in different spectral states, as monitored by the Rossi X- 
ray Timing Explorer over about a month. We fit the power density spectra 
using multiple Lorentzians. We show that the characteristic frequencies of these 
^ ■ Lorentzians, when properly identified, fit within the correlations previously re- 

Q"^ ■ ported. The time evolution of these frequencies and their relation with the pa- 

^ ■ rameters of the energy spectra reported in Barret & Olive (2002) are used to 

O ■ constrain the accretion geometry changes. The spectral data were fitted by the 

^ ■ sum of a blackbody and a Comptonized component and were interpreted in the 

framework of a truncated accretion disk geometry, with a varying truncation ra- 
dius. If one assumes that the characteristic frequencies of the Lorentzians are 
■ some measure of this truncation radius, as in most theoretical models, then the 

timing data presented here strengthen the above interpretation. The soft to hard 
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and hard to soft transitions are clearly associated with the disk receding from and 
approaching the neutron star respectively. During the transitions, correlations 
are found between the Lorentzian frequencies and the flux and temperature of 
the blackbody, which is thus likely to be coming from the disk. On the other 
hand, in the hard state, the characteristic Lorentzians frequencies which are at 
the lowest, remained nearly constant despite significant evolution of the spectra 
parameters. The disk no longer contributes to the X-ray emission, and the black- 
body is now likely to be emitted by the neutron star surface which is providing 
the seed photons for the Comptonization. 
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1. Introduction 

Rapid X-ray variability is a powerful probe of physics of accretion flows around neutron 
stars and black holes. Significant progress has been accomplished recently with the advent 
of the Rossi X-ray Timing Explorer [RXTE, Bradt, Rothschild & Swank 1993). Its large 
collecting area and micro-second time resolution allowed us to study with exquisite details 
the rapid variability of a wide variety of accreting X-ray sources. In particular, the power 
density spectra (PDS) of Galactic black hole candidates and neutron star low-mass X-ray 
binaries (LMXBs) exhibit a variety of features ranging from narrow quasi-periodic oscillations 
(QPOs) to broad noise components (for a review see Wijnands 2001 and references therein). 

It has been shown recently that these PDS can be well represented as a superposition of 
a few Lorentzians (e.g. Olive et al. 1998; Van Straaten et al. 2002; Belloni, Psaltis & Van der 
Klis 2002). The multi-Lorentzian approach gives a simple and universal phenomenological 
description of the PDS. Each noise component (Lorentzian) is then described by its charac- 
teristic frequency, width and amplitude, usually expressed as a fractional root-mean-square 
(RMS). It has been demonstrated that these characteristic frequencies are correlated for a 
given source and between sources, both for neutron star and black hole binaries (Wijnands 
& Van der Klis 1999; Psaltis, Belloni & van der Khs 1999; Belloni et al. 2002). In most 
models of rapid X-ray variability, these frequencies are related to some kind of an abrupt 
transition in the accretion disk (e.g. Stella, Vietri & Morsink 1999; Titarchuk & Osherovitch 
1999; Psaltis & Norman 2002). This could be a truncation of a Shakura-Sunyaev disk into 
a hot inner flow (e.g. Rozahska & Czerny 2000). With the truncation radius decreasing (the 
disk coming in) the characteristic frequencies increase. 

4U 1705-44 is a neutron star LMXB system which belongs to the class of atoll sources 
(Hasinger & Van der Klis 1989). The source shows variability on all time scales, from months 
down to milliseconds (Langmeier et al. 1987; Berger & Van der Klis 1998, Ford, van der Klis 
& Kaaret 1998, Liu, van Paradijs & van den Heuvel 2001 and references therein). On long 
time scales, it displays clear luminosity-related spectral changes. This is illustrated in the 
observations reported by Barret & Olive (2002; hereafter Paper I) who followed the spectral 
evolution of the source during a transition between soft and hard X-ray spectral states. The 
data were interpreted in the framework of a truncated accretion disk geometry for which the 
spectral transitions observed were associated with changes in the disk truncation radius. 
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In this paper, we follow the multi-Lorentzian approach to study the rapid X-ray vari- 
ability of the source using the same set of observations as presented in Paper I. The aim 
of this work is to look for correlations between the timing and spectral parameters to get 
further insights on the accretion flow changes associated with the state transitions. 



2. Observation and data analysis 

We use the same set of RXTE observations as in Paper I. The data cover a period 
from February 10th to March 9th, 1999, with pointings performed every two days, giving 
14 observations in total. In the following the observations are numbered from 1 to 14. 
During this period, the source X-ray count rate first decreased, reached a minimum around 
observation 6, and then increased to reach the initial level (sec top of Fig. 1). On the X-ray 
color-color diagram (bottom of Fig. 1), the source traced out part of a more or less triangular 
pattern (see also Gierlinski & Done 2002). 4U 1705-44 evolved from the soft (the so-called 
banana) to the hard (island) state in ~ 10 days, remained in a hard state for ~ 12 days, and 
then moved back to the soft state in ~ 4 days. In Paper I, all fourteen energy spectra were 
adequately fitted by the sum of a blackbody, thermal Comptonization and an iron line. 

For timing analysis we use the science event data obtained in configuration E_16us_64M_0_ls. 
We have computed 0.03-4096 Hz PDS in the 3-30 keV energy range, for the same time in- 
tervals as the energy spectra. The Poisson counting noise level was estimated from the data 
above 1500 Hz and subtracted from the PDS before fitting. Each PDS has been rebinned 
using a logarithmic scheme, and expressed in terms of fractional RMS amplitude. 

We have fitted each of the 14 PDS by up to 4 Lorentzians with the following analytical 
formula given below (van Straaten et al. 2002): 



(f + arctan 2Q) [A2 + (z/ - 2^Qf] ' ^ ' 

where A = fmax/\/ 1 + 4Q2 is the half width at half maximum of the Lorentzian and i/max (the 
position of the maximum of the Lorentzian per logarithmic frequency interval, i.e. the peak 
position of vPv)^ r (the fractional RMS integrated from to oo) and Q (the quality factor) 
are the independent parameters to fit. The ccntroid of the Lorentzian (z/q) is related to the 
above parameters as: vq = 2AQ =\/^^max — A^. The spectral fitting and error computations 
were carried out using XSPEC 11.1 (Arnaud et al. 1996). Whenever the Q parameter was 
consistent with zero, we used zero-centered Lorentzians {Q frozen at 0). Finally, we searched 
for kHz QPOs using the PDS scanning technique described in Boirin et al. (2000). 
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3. Results 

The spectra together with the best-fitting models are presented in Fig. 2. The best fit 
parameters are hsted in Tab. 1. Time evolution of the fitted Lorentzian frequencies is shown 
in Fig. 3. In those observations with the best statistics (7 to 13), three broad components are 
detected with noise extending up to > 500 Hz. kHz QPOs around 750-800 Hz were detected 
above the 3a level in observations 1, 2 and 14 (see Tab. 2) when the source was in the soft 
state. Our frequencies are close to those reported by Ford et al. (1998) using observations 
performed in February to June 1997 while the source was in a similar spectral state. In 
terms of fractional RMS, the sensitivity to QPO detection decreases with the source count 
rate. For indication, Sa upper limits of 4.6,8.5,6.2,5.0,4.5,3.4% can be derived for a QPO 
at 750 Hz with a 20Hz width in observations 3, 5, 7, 9, 11, 13 respectively. However, in 
terms of modulated amplitude, the signals detected in observations 1, 2 and 14 would have 
been clearly detected in the other observations. 

Narrow QPOs around ~ 170 Hz were found in observations 2 and 3 {u = 162^^ Hz, 
RMS(%)=2.9 ±0.7 for observation 2 and z/ = UStlo Hz, RMS(%)=4.5 ± 1.0 for observation 
3). Similar QPOs, called hecto-Hz QPOs have been reported from other sources (e.g. van 
Straaten et al. 2002). The PDS of observation 2 is also remarkable by the presence of a very 
low frequency component (VLFN, see Fig. 2), not seen in other observations. A fit with 
a zero-centered Lorentzian of this component yields u^ax — 0.156lo;o4 Hz and RMS(%) = 
3.1 ± 0.4. These transient features will not be discussed in this paper. 

3.1. Identification of the PDS components 

The PDS of 4U 1705-44 resemble those reported by Belloni et al. (2002) and van Straaten 
et al. (2002) from similar systems. Following Belloni et al. (2002), we name the first two 
Lorentzians as band limited noise (BLN at uj, ~ 1-7 Hz) and low-frequency noise (LFN 
at i/LF ~ 7-40 Hz). The LFN component is the only one which can be identified in all 
observations. In Fig. 4, we plot z/lf versus z/5. A clear correlation can be seen. Our results 
can be directly compared with those obtained for some other neutron star LMXBs. Taking 
out the same two components in Belloni et al. (2002) for IE 1724-3045 and GS 1826-34 and 
van Straaten et al. (2002) for 4U 1728-34 and 4U 0614-^091 (in the latter paper Ub is called 
i^bln), we find that our data points fall in on the global correlation previously reported by 
these authors. This supports the identification of the first two components of the PDS of 
4U 1705-44 with BLN and LFN. 

The third broad Lorentzian component which we denote as high frequency noise (HFN 
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at characteristic frequencies of ~ 60-400 Hz) is detected only in data sets 7 to 13. It is the 
one which is the less statistically significant, and its identification is therefore less certain. 
It may be the hecto-Hz component reported by van Straaten et al. (2002) from similar 
systems. However, these hecto-Hz components arc generally narrower than the features we 
observed in 4U 1705-44 which are fitted with zero-centered Lorentzians. It seems therefore 
more likely to be what is called L/ in Belloni et al. (2002); namely the extrapolation towards 
low frequencies of the lower kHz QPO. Support from that hypothesis comes from the fact 
our data points match the global correlation between i/lf and i/i discovered by Psaltis et 
al. (1999). This feature would then be closely related to the lower kHz QPOs detected in 
observations 1, 2 and 14. As expected, the kHz QPO frequencies fit in the global correlation 
(see Fig. 5). 



3.2. Correlation between PDS and energy spectrum parameters 

Since the LFN component is detected in all observations and since the other frequencies 
are correlated with i/lf, we used i^lf for searching correlations between PDS and energy 
spectrum parameters. 

We split the observation in 3 segments: the transition the soft to hard state (Tl, obs. 
1 to 6), the observations corresponding to the island state (IS, obs. 6 to 12), the transition 
hard to soft (T2, obs. 12 to 14). In Fig. 6, 7 and 8, we plot i/lf versus the soft and hard 
colors, i^LF versus the blackbody temperature and bolometric fiux and z/lf versus the electron 
temperature of the Comptonizing cloud (kTe) and 20-200 keV hard X-ray fiux (Fnard)- 

Fig. 6 shows that there is a clear anti-correlation between j/lf and HR2 in all three 
segments. On the other hand, there are no obvious correlations between z/lf and HRl. This 
is consistent with previous work which have shown that timing features, such as kHz QPOs 
(whose frequencies correlate with uj^p, see above) anti-correlate with the hard color (e.g. 
Mendez et al. 1999). 

Fig. 7 clearly shows a nice correlation between l/lf and the blackbody temperature and 

flux at the start of the soft to hard transition (obs. 1 to 3) and during the reverse transition 
(obs. 12 to 14). Similar correlations between the blackbody flux and kilo-QPO frequencies 
have been already reported from 4U 0614+091 (Ford et al. 1997). Here we show that the 
correlations break down in 4U 1705-44 when the blackbody flux dropped below ~ lO^^'' ergs 
s~^ cm~^. Down in the IS, the absence of correlation is striking, despite the smooth but 
significant increase of the blackbody fiux and temperature. 

Finally, Fig. 8 shows that during Tl and T2 there is a clear anti-correlation between 
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i^LF and both kTe and Fnard, and there are no obvious correlations of these two parameters 
with i^LF in IS (no correlations with time are present either). 

4. Discussion 

In most of the rapid variability models the frequencies of particular noise components 
are set at some transition radius in the disk (e.g. Stella et al. 1999; Psaltis & Norman 2002). 
It could be for instance a truncation of a standard Shakura-Sunyaev disk into a hot inner 
flow. Here we assume that the frequencies wc observe, in particular z/lf which can be traced 
throughout all the data sets, are related to the truncation radius of the disk. In Paper I 
we interpreted the spectral state transitions in the framework of a such truncated accretion 
disk geometry, allowing the truncation radius to change within the observations. Let us now 
discuss the new constraints that the data presented above set on this scenario. 

At the beginning of our observations (obs. 1-2), the source is in the soft state, the 
truncation radius is small and the disk is close to the neutron star surface (see Paper 1). The 
characteristic frequency z/lf is at its maximum around 40 Hz, which is consistent with the 
inner disc radius reaching its minimum. It is also the time during which the kHz QPOs at 
750-780 Hz are detected. If the latter one is associated with the Keplerian frequency at the 
truncation radius, the disk truncates at ~ 20 km. The soft component (/cTrb ~ 1-8 keV) 
present in the energy spectrum originates from the disk (Paper 1). The hard component is 
from the boundary layer, with low temperature of kT^ ~ 4 keV. 

In Paper I, we suggested that during the transition to the hard state (Tl, obs. 2 to 6) 
the disk receded while the total source luminosity decreased from ~ 2.1 to 0.7 x lO^^ergs 
s~^. This is supported by our timing analysis: during the transition i/lf dropped roughly 
by factor of ~ 6. The fastest variability timescale around a compact object is given by the 
Keplerian orbital frequency (at v^). This sets an upper bound on the inner radius of the 
disk at which the variability is produced (see e.g. di Matteo & Psaltis 1999). When is 
not detected, one can use the empirical correlation found by Psaltis ct al. (1999) between 
and i^LF to estimate the inner disk radius as a function of z/lf (see equation 2 in di Matteo & 
Psaltis 1999). Using the latter equation, one obtains a disk truncation radius around 70 km 
for observation 6. When the disk recedes, the temperature and flux of the soft component 
drop to their lowest values (kTBB~ 0.8 keV, Fig. 7, left panels). Meanwhile, the hot inner 
flow (or optically thin boundary layer) builds up and its temperature rises from ~ 4 to ~ 14 
keV (Fig. 8, left panels) possibly because the disk contribution to the cooling of the inner 
flow is then seriously reduced. 
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Between observations 6 and 12 (IS) the source is in the hard state while the mass accre- 
tion rate (inferred from the source luminosity Lo.i_2oo ~ 0.7 to 2.5 x lO^^ergs s~^) monoton- 
ically increased. The energy spectrum is dominated by the emission from the hot optically 
thin inner flow (Paper 1). In the proposed scenario, the disk is far away, and z/lf is low, ~ 
7-12 Hz. With the exception of observation 6, all the PDS show three broad Lorentzians (at 
Ub, z/LF and ui) at characteristic frequencies which are roughly constant suggesting that the 
geometry of accretion remains similar through these observations. There is no correlation 
between i/i,p (which is roughly stable) and the blackbody flux and temperature either (Fig. 
7, right panels) despite significant variations of these parameters. During this hard state, 
the spectral and the timing parameters are de- correlated. This can be naturally explained 
if the soft component (blackbody) originates from the neutron star surface, not from the 
disk. Basically, the disk is far away and its temperature is so low that it cannot be seen 
by the PCA, so it does not contribute to the energy spectrum. We can see its variations 
only through the timing properties, as z/lf- The blackbody comes from the neutron star sur- 
face which is continuously heated along these observations. This soft component is further 
Comptonized in the optically thin hot inner flow/boundary layer leading to an increasing 
-^20-200 flux (Fig. 8, right panel). It is also important to note that at the end of the IS state 
(obs. 12), the source has recovered an accretion rate (as inferred with the luminosity) similar 
to observations 1 or 2, while the truncation radius of the disk (as traced with i/lf) is still 
large. 

The transition back to the soft state (T2, obs. 12 to 14) is the inverse of the Tl transition. 
The blackbody flux and temperature increased, which was interpreted (Paper I) as the 
reappearance in the PCA energy band of the disc moving in and becoming hotter. Again, 
the timing analysis supports this result: there is a significant increase in z/lf, very well 
correlated both with Tbb and Fbb- Meanwhile, the inner flow is significantly cooled by the 
increasing flux of soft photons coming from the disk and its temperature falls from ~ 10 to 
~ 4 keV. In observation 14, the system is back to a state very similar to observation 1 and 
2, and we start seeing again kHz QPO at ~ 750 Hz. 

It is also interesting to look at rapid X-ray variability as a function of the position in 
the color-color diagram. Both atolls and Z sources trace out a characteristic track in the 
diagram (see e.g. Hasinger & van der Khs 1989). The frequencies in the PDS are usually 
very well correlated with the position in the diagram. It has been shown recently that at low 
luminosities atolls form an additional, upper horizontal branch in the color-color diagram, 
tracing out a Z-shaped track, similar to Z sources (Gierhnski & Done 2002; Muno, Remillard 
& Chakrabarty 2002). 

4U 1705-44 moved from left to right in the horizontal branch between observations 6 and 



-8- 



12, and then along the diagonal, between observations 12 and 14 (Fig. 1). The timing and 
spectral properties on the diagonal branch are typical of other atoll sources. However, the 
upper branch is different: as we discussed previously, there is no correlation between z/lf and 
the position on the branch, which in this case (branch is horizontal) can be described by the 
soft color HRl (see Fig. 6). 

In our scenario, at low accretion rate, the position in the upper branch of the diagram 
is set by the neutron star temperature, which also supplies the seed photons. With the 
increasing accretion rate, M, this temperature increases, so does the temperature of the 
low-cncrgy cutoff in the Comptonizcd component. The spectrum below 4 kcV hardens, so 
the soft color increases and the source moves to the right in the diagram. This movement is 
not accompanied by the decrease of the truncation radius in the disk. Then, on the diagonal 
branch the disk temperature increases, so it can be seen in the PCA spectrum, while the 
neutron star surface is obscured by the increasing optical depth of the boundary layer. Thus, 
both spectral and timing properties depend now on the disk, and therefore are correlated. 
In this picture, the position on the upper branch depends on M, while the characteristic 
frequency i/lf do not. 

The path of 4U 1705-44 in the color-color diagram is even more comphcated. On 

its way from the banana to the island state (Tl) it traces out a different track than in the 
opposite direction (T2). The timing and spectral properties on both tracks arc different. For 
instance observations 3 and 13 have very similar uj^p, but different bolometric luminosities 
(Lo.i-200 = 1-4 X 10'^'' ergs s~'^, 2.4 x 10^"^ ergs s^^ respectively, sec Paper 1) and different power 
density spectra (see Fig. 2). This shows that for the same truncation radius, one can have 
different instantaneous accretion rate as derived from the source luminosity, and different 
time variability (the Lorentzians at and are not detected in observation 3). In the 
proposed scenario, there are however important differences between observations 3 and 13 
(and more generally between Tl and T2). In the first case, the disk recedes, the luminosity 
is decreasing, whereas in obs. 13 the disk approaches the neutron star and the luminosity 
is increasing. This means that the state of the system depends at some level on its recent 
history. Any models proposed to explain the changes in z/lf must take this fact into account. 

One possibility could be that z/lf (i-e the truncation radius) is set by some long-time 
averaging process over time scales of days and only depends on a parameter 77 = M^/ <Md> 
(where Md is the accretion rate through the disk) similarly to what has been proposed by van 
der Klis (2001) to explain the parallel tracks phenomenon in LMXBs. This model naturally 
predicts the decorrelation between the truncation radius position and the total accretion rate 
(and so the luminosity). In that scheme, the motion of the disk is associated with a rapid 
decrease of Md during the Tl transition (i.e 77 < 1, decreasing with time, the disk recedes) 
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and a reversed evolution during T2. The disk is stable either if Ma is constant or if it varies 
linearly with time over several days. This could be the case from observations 6 to 12 when 
i^LF was low and roughly stable. 



5. Conclusions 



The spectral and timing data of the state transitions observed from 4U 1705-44 can be 
explained in the framework of a model, in which the critical parameter is the position of 
the truncation radius between the disk and a hotter inner flow. This parameter sets both 
the frequencies of the timing features and the spectral shape. What sets the value of the 
truncation radius is unclear at the moment, but both the timing and spectral data indicate 
that it cannot be the mass accretion rate derived from the bolometric source luminosity. One 
possibility could be this radius is set by some long-time averaging process over time scales 
of days. 
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Table 1: Best fit parameters for the 14 PDS. z/^ and RMSb are respectively the frequency 
and the RMS of the band-limited noise component fitted with a zero-centered lorentzian. 
i^LF RMSlf and Qlf are the frequency, the RMS and the Quality factor of the Low- Frequency 
Noise component detected in all observations, ui and RMSi are the frequency and RMS of 
the high frequency noise componant. The naming convention follows Belloni et al. (2002). 
Finally, /d.o.f is the value and the number of degree of freedom of the fit. 
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Table 2: The best fit parameters (frequency, RMS and the quality factor Q) for the kHz 
QPOs detected in observations 1, 2 and 14. For observation 14, the Q value was poorly 
constrained and then frozen at 10 for the fit. These frequencies are consistent with those 
reported by Ford et al. (1998). 
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Fig. 1. — The PCA background subtracted light curve of the 4U 1705-44 in the 3-16 keV 
band (top panel). The data are plotted versus a time index to avoid the gaps in time between 
observations. The bottom panel shows the X-ray color-color diagram of 4U 1705-44. The 
soft color (HRl) is defined as the ratio between the 4.3-6.5 keV counts and 2.9-4.3 keV 
counts, and the hard color (HR2) as the ratio between the 10.1-16.2 keV counts and 6.5-10.1 
keV counts. The observation numbers allow to follow the path of the source during the state 
transitions. 
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Fig. 2. — The 14 power density spectra (PDS) from 4U 1705-44 fitted by up to four 
Lorentzians. HFN stands for High Frequency Noise, LFN for Low-Frequency Noise, BLN 
for Band-Limited Noise and VLFN for Very-Low Frequency Noise. Narrow Quasi-Periodic 
OsciUations (QPOs) were also detected in observations 1, 2, 3 and 14. 
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Fig. 3. — Top panel: The luminosities {LTotai, Luard and LBiackhody) along the 14 observations 
in units of lO^^ergs s~^. Bottom panel: Time evolution of the PDS fitted parameters, z/j is 
the frequency of the lorentzian fitting the BLN component of Fig. 2, z/lf is the frequency of 
the LF lorentzian, vi is the HFN Lorentzian. The naming convention follows Belloni et al. 
(2002). Hecto-Hz QPOs were detected in observations 2 and 3. kHz QPOs were detected in 
observations 1, 2 and 14. As can be seen i/lf is the only component which is detected in all 
observations. 
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Fig. 4. — i/LF as a function of z/^. Data points from 4U 1705-44 are plotted with a filled 
circle. They are compared with those obtained for systems similar to 4U 1705-44 for which 
the same type of multi-lorentzian fitting of the PDS was performed (Empty stars : GS 1826- 
34, IE 1724-3045 and SIX 1735-369; Belloni et al. (2002). Empty crosses : 4U 0614+091 
and 4U 1728-34; van Straaten et al. (2002)) Our data point fit in the correlation, indicating 
that the identification of the various frequencies of the PDS is correct. 
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Fig. 5. — i^LF as a function of vi. For comparison data points from Belloni et al. (2002) are 
also shown. Again, a correlation has been shown to exist in many systems (Psaltis et al. 
1999; Belloni et al. 2002). vi and the three kHz QPO frequencies follow the correlation. 
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Fig. 6. — i/LF as a function of HRl and HR2. Three segments of the observations are 
considered: Tl is the transition which occurred between observations 1 and 6 (filled circles), 
IS is the Island State between observations 6 and 12 (empty squares) and T2 is the transition 
between observations 12 and 14 (empty circles). 
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Fig. 7. — i^LF as a function of the blackbody temperature and flux for the three segments of 
observations. 
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Fig. 8. — i^LF as a function of the electron temperature and the bolometric flux between 20 
and 200 keV for the three segments of observations. 



